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Interference of Light in Michelson— Morley Interferometer: 
A Quantum Optical Approach 
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We investigate how the temporal coherence interference properties of light in a Michelson-Morley 
interferometer (MMI), using only a single-photon detector, can be understood in a quantum-optics 
framework in a straightforward and pedagogical manner. For this purpose we make use of elementary 
quantum field theory and Glaubers theory for photon detection in order to calculate the expected 
interference pattern in the MMI. If a thermal reference source is used in the MMI local oscillator 
port in combination with a thermal source in the signal port, the interference pattern revealed by 
such an intensity measurement shows a distinctive dependence on the differences in the temperature 
of the two sources. The MMI can therefore be used in order to perform temperature measurements. 
A related method was actually used to carry out high precision measurements of the cosmic micro- 
wave background radiation on board of the COBE satellite. The theoretical framework allows us 
to consider any initial quantum state. The interference of single photons as a tool to determine 
the peak angular-frequency of a one-photon pulse interfering with a single-photon reference pulse 
is, e.g., considered. A similar consideration for laser pulses, in terms of coherent states, leads to a 
different response in the detector. The MMI experimental setup is therefore in a sense an example 
of an optical device where one can exhibit the difference between classical and quantum-mechanical 
light using only intensity measurements. 
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I. INTRODUCTION 

In 2006 G. Smooth and J. Mather shared the No- 
bel Prize in Physics "for their discovery of the black- 
body form and anisotropy of the cosmic micro- wave back- 
ground radiation (CMB)" These exciting discoveries 
were a breakthrough in modern cosmology by the CMB 
anisotropy and the strong validation of the black body 
spectrum as predicted by the Big Bang theory. The dis- 
covery of the black body form of the CMB spectrum and 
the high precision measurement of the CMB tempera- 
ture (see e.g. Ref.Q ) relied heavily on the so called 
Far Infrared Absolute Spectrophotometer (FIRAS) H on 
board the Cosmic Background Explorer (COBE) [J,Q- 
In short the FIRAS is a Michelson-Morley interferom- 
eter enabling a comparison of the interference patterns 
between an observed source and a reference source on- 
board the COBE satellite. 

In this paper we will make use of Glaubers theory for 
photon detection [1,01 (for a guide to the early liter- 
ature see e.g. Refs.Q and for text-book accounts see 
e.g. Refs.|9l4ll|) together with elementary quantum me- 
chanics to show how the principles of the FIRAS can be 
understood in a quantum-optics framework. 

In Section we |TT] recapitulate the principles of Glaubers 



photon detection theory and the transformation prop- 
erties of a quantum field in a beam-splitter (see e.g. 
Refs. [l2l - [l3 |). The Glauber theory of optical coher- 
ence is by now well established and plays a central role 
in fundamental studies of quantum interference effects 
of photon quantum states (see e.g. Refs. [l5l - fr7j ). In 
Section IIIII we consider temporal interference effects in 
the Michelson-Morley interferometer for pure quantum 
states like single-photon states as well as classical states 
corresponding to coherent states (see e.g. Refs.Q). In 
Section IIVI we explain the principles of interference of 
thermal light in the Michelson-Morley interferometer by 
using only vacuum as the reference source and reproduce 
known expressions and, in Section [Vj we consider the full 
system with an observed thermal source combined with 
a thermal reference source. With the results obtained we 
then explain the basic principle of FIRAS and how it was 
used as a high precision thermometer. In Section fVTl we. 
finally, give some concluding remarks. 

Our presentation extends a recent presentation by 
Donges [l8| and illustrates, e.g., that a quantum- 
mechanical treatment directly leads to the concept of a 
thermal coherence length without explicitly making use 
of classical results like the Wiener-Kintchine theorem as 
in Ref.[l||. 
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II. THE MICHELSON-MORLEY 
INTERFEROMETER 

We consider the Michelson-Morley interferometer 
(MMI) as illustrated in FigJTJ where the so called tempo- 
ral coherence properties of the radiation field is probed 
(for an early account see e.g. Ref. 19]). In order to under- 
stand the appearance of interference effects in the MMI, 
we first discuss the separate parts of the MMI before we 
consider the full setup with the presence of a reference 
beam. 



A. Glaubers theory of photon detection 

Let us first outline a simple, but not unrealistic, model 
of a photon detector situated at the space-time point 
(x, t) . In this simplified model of a photon detection pro- 
cess [1, 0] , the detection of the photon is described by an 
annihilation of a photon at the detector which modifies 
the initial state as follows: 



|t) E i+) (x,t)\i) 



(1) 



The observable electric field operator E(x,t) = 
^j m E m {x,t) is described in terms of a suitable normal 
mode expansion, indexed by mode the number m, as a 
superposition of positive and negative angular-frequency 
contributions: 



E m (x,t) = Eit\x,t) + E^\x,t) , 



(2) 



where Em\x,t) (Em\x,t)) contains an annihilation 
(creation) operator for a photon with mode number m. 
According to the basic Born rule in quantum mechanics, 
the probability to detect the system in a final state |/), 
after the one-photon absorption process, is then propor- 
tional to |(/|i?(+)(£,£)|i)| 2 . Since the exact details of the 
final states are, in general, unknown we sum over all pos- 
sible final states |/). In general, we also have to consider 
not only a pure initial quantum state but also a quantum 
state as described by a density matrix. This leads to a 
description of the observed intensity / which we write in 
form 



I = Tr IpE^ (x, t)E {+) (x, t) 



(3) 



where p is the density matrix describing the initial state, 
and where use have be made of the completeness of all 
possible final states, i.e. J2f \f)(f\ = 1- Since we will 
not be interested in the absolute normalization of the 
observed intensity J, we can neglect normalization con- 
stants that may enter into /. It is a remarkable fact that 
an analysis of single-photon interference in a Young in- 
terferometer using such a quantum-mechanical descrip- 
tion of the photon detection process is fairly recent in 
the history physics [20| as well as a proper experimental 
investigation of single-photon interference [U [22j ■ 



For reasons of clarity, we will now consider a normal 
mode expansion of the electro-magnetic field observable 
E(x, t) in terms of plane waves, i.e. in terms of a mode 
sum over wave- vectors k and polarization- vector labels A, 
i.e. we have for its positive angular-frequency part that 



E [+) (x, t) = i 



E 



ik-x — iuifct 



(4) 



Here w& = c\k\ and t is a suitable retarded time param- 
eter which will be given in terms of the time-of-flight, 
given the optical paths in the MMI setup. If the direc- 
tion of the light beams considered are well-defined, the 
dependence of the detector position x can be neglected 
in the expression for I. In general spatial modulations 
of the measured intensity are expected [23|, [24| . A the- 
oretical analysis of such effects alon g lin es as discussed 
in the literature (see e.g. Refs. (25, 26]) will, however, 
not enlighten the issues we are addressing in the present 
paper. V is a quantization volume that will be allowed 
to be arbitrarily large at an appropriate late stage of our 
calculations. is the annihilation operator describing 
the detected light and eV, denotes the unit polarization 
vector of the normal mode considered. 

Since, in the end, the dependence of normalization con- 
stants will be irrelevant, and since we will only consider 
polarization-independent optical devices, we make use 
of a scalar notation. We therefore suppress the wave- 
vector and polarization labels and with u) = uik we write 
a(cj) = a£ A , such that [a(u>), <v(cj)] = S uu i in terms of 
a discrete Kronecker delta 8 UU >. We also make use of 
the following convenient notation [n| for the positive 
angular-frequency part of the electric field at the posi- 
tion of the detector at time t 



£<+>(*) = iJ2 ^ 1/2 «(^)e #(i) 



(5) 



Here <j>(i) = — ui(t — r s ) now denotes an optical phase 
which explicitly takes the source-detection retardation 
time into account by the time-delay r s , which will be 
evaluated for the MMI-setup below. If the detector time 
t enters explicitly into the detection intensity Eq. (J3j , we 
will perform a time-average which corresponds to a fi- 
nite detector time-resolution window. The correspond- 
ing time-average of the observed intensity / = I(t) will 
be denoted by (I), i.e. 



(I) 



T 



7W2 



dtl(t) . 



(6) 



int J -T int / 2 



where the time Ti„ t of integration, as e.g. the time during 
which an actual measurement proceeds, is chosen to be 
sufficiently large in comparison with any reciprocal band- 
width of the initial quantum states |z) considered. The 
time-averaged observed intensity (/) will in general, as 
we will see explicitly below, be a function of a time-delay 
r depending on the actual experimental set-up. 
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For a finite quantization volume V, ui can be regarded 
to be discrete and, in the infinite volume limit, we make 
use of the rule that XL ~~ > Jo° duiui d ~ 1 , where d is the 
number of space-dimensions. We will neglect the depen- 
dence of transverse dimensions of very collimated normal 
modes and assume that d = 1. When appropriate we will, 
however, also consider d = 3 in order to compare with 
related results in the literature [18j, ll9j . Our main results 
will, however, not be very sensitive to the choice of d. 



Detector 




Mirror 1 

FIG. 1: Schematic drawing of the MMI setup considered in 
this paper. The thick and diagonal line is representing a 
beam-splitter with a transmission coefficient T. The numbers 
in the parenthesis are referring to the indices use to describe 
the different light beams in the calculations. The signal port 
(1) and the local oscillator port LO (0) are prepared with 
various quantum states, like Fock states, coherent states or 
states with random phases like thermal states. In order to 
describe photon absorption in the detector we make use of 
Glaubers photon detection theory. 



Using the same notation as above, a one-photon quan- 
tum state |/), with an angular-frequency distribution 
given by / = /(w), is then given by 

|/) = ^^/HIU^(/, a t)|o), (7) 

where |l w ) = a^w)^) denotes a one-photon state with 
angular-frequency uj, normalized according to (l w |lcy) = 
6 uu i, and where |0) denotes the vacuum state. We also 
make use of the notation (/, aJ) = J2 U Vfiw f (w)aJ (m) . 
The state |/) above is an eigenstate of the number opera- 
tor N = J2uj o,^ (uj)a(uj), i.e. a Fock state, with, of course, 



an eigenvalue corresponding to one particle present. Nor- 
malization of the state |/) for d = 1 therefore corresponds 
to 

(/|/) = i = E<H/HI 2 = / cM/Ml 2 , (8) 

in the large- volume V limit. In order to make our pre- 
sentation quantitative we will, for reasons of simplicity, 
consider real- valued one-photon angular-frequency distri- 
butions /(w) such that 

/H = lexp(-( W -^) 2 /2 ( r 2 ), (9) 

with a mean angular-frequency ui and width a and where 
the normalization constant N is given by 

'"''^(^f (10) 

in terms of an error function. This choice of frequency 
distribution makes it possible to actually carry out all 
relevant expressions analytically. In obtaining the prop- 
erly normalized expression Eq.® we keep uj > 0. It 
may, however, sometimes be possible to extend the range 
of angular frequencies to arbitrarily negative values in 
Eq.®, so that |iV| 2 = air, with an exponential small 
error, which makes some of the expressions as given be- 
low more tractable and transparent. With our choice 
of beam parameters below it turns out that such an ap- 
proximation plays only a minor role with regard to actual 
numerical evaluations. 

Conventional coherent states |/) c , as expressed in 
terms of the one-photon distribution /, can then be ob- 
tained using a multi-mode displacement operator (see e.g. 
Refs.ll), Le - 

|/) c = e (/^)-(r,a)| 0) = e -(/|/)/2 e (/,at)| 0) _ (n) 

such that 

aM) c = (5") 1/2 m\f) c , (12) 

and hence 

c (/i^i/> c = Eh/hi 2 - (13) 



B. Transformation in the beam-splitter 

Next, we consider a beam-splitter with frequency inde- 
pendent transmittance T and reflectance R. If we assume 
a prefect beam-splitter, where all light is cither reflected 
or transmitted, we have R + T = 1. The input annihila- 
tion operators ao(w) and a\(u>) of the bea m-sp litter will 
then transform according to (see e.g. Refs. |12Hl4j ) 

02(0;) = VTao(u>) + iy/T^-Tai(uj) , 

a 3 (uj) = Vr ai (uj)+iVl-Tao(uj), (14) 
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where oq(uj), ai (to) are the L.O. and signal port mode an- 
nihilation operators, and 02(10), 0,3(01) the output annihi- 
lation mode operators corresponding to the transmitted 
and reflected modes, respectively. The 7r/2-phase-shift 
between the transmitted and reflected part, described by 
the complex numbers in Eq. (|14|) will play and important 
role below as is also the case in the famous Hong-Ou- 
Mandel two-photon correlation experiment |27| and re- 
lated investigations (see e.g. Ref. [2la - [3fj| | ) . 



1.10 




FIG. 2: The normalized single-photon intensity (7)(t)/(7)(0) 
as a function of the dimensionless time-delay to for the case 
of one-photon states in the signal and the LO ports with the 
same spectral width a but with different mean frequencies. 
Two different examples are plotted with ui a = 3<r: uji — 3.15a 
with (solid line) or uii = 2.85<r (dashed line). The asymptotic 
values of (7) (r) / (7) (0) can be obtained from the expression 
Eq. (|2ip in the main text, i.e. (1 + uii /ui 3 )/2. 



A light beam arriving at the beam-splitter after being 
reflected in the mirrors 1 and 2 will, in general, be phase 
shifted, i.e. expressed in terms of mode operators this 
process corresponds to the time-of-flight replacement 



a 2 M 

03 M 



a 2 (u;)e l<fc(t) 
a 3 (u;)e l<fe(t) 



(15) 



due to difference in optical path lengths with <f>2 — 
—ui(t — T2), 4>3 = — ui(t — T3) in terms of time-delays 
T2 and T3. The reflections at the mirrors in the MMI 
setup will also introduce a phase-shift of 7r/2, but this is 
equal for the two light beams and we can therefore be 
neglected all together. A light beam passing through the 
beam-splitter after reflection at the mirrors will, again, 
transform according to Eq. (fl"4|) and we therefore, finally, 
obtain an expression for the mode operator describing 
incident light on the photon detector, i.e. 

a 5 (co) = VTa 2 (w)e* (t) + iVT~^ Ta 3 (w)e^ 3(t) 
= ao(w) (Ye^ 2(t) - (1 - ly^w) + 



oi(«) (iy/T(l - T)e*« + iy/T(l - T)e 



# 3 (t) 



(16) 



The positive angular-frequency part of the electro- 
magnetic field operator, E^(t), at the position of the de- 
tector as given in Eq. ([5]) , will now be expressed in terms 
of 05(0;). For reasons of clarity we will consider a 50/50 
beam splitter, i.e. we make the choice T = 1/2, and 
therefore the corresponding electro-magnetic field opera- 
tor to be used in Glaubers theory of photon detection, is 
given by 



£(+)(*) = ^Vd^ui 1 / 2 (aoM^" - e*»«)+ 



Wi(w)(' 



,itp2(t) 1 e i<t>3(t) 



(17) 



III. 



INTERFERENCE OF FOCK STATES AND 
COHERENT STATES 



Let us now specifically consider the following initial 
Fock state 



i) = \fs) ® \fio) , 



(18) 



for the signal and the local port with one-photon angular- 
frequency distributions according to Eq. © with the same 
spectral widths a but with Co = ui s = 3<r for the signal 
port and Q = uj; = 3.15cr or Cdi — 2.85ct for the local 
oscillator port. We then see that 

E {+) m) = 

i £ ^Slooj 1 / 2 (i/ s ( W )e*W(l + e^)|0) ® |/ /o ) 

UJ 

+/ io ( W ) e *( t )(i- e — ))\f s )®n) , (19) 



where we have made use of the fact that 4>$ (t) — cf>2 (t) = 
lo(t3 —T2) = lot, which now defines the optical time-delay 
t. The probability for the detection of one photon will 
according to Glaubers theory of photon detection, as we 
have mentioned above, be proportional to the absolute 
modulus square of the probability amplitude E^ + \t)\i), 
i.e. to |i?( + )(i)|?)| 2 . A time-average over the time t ac- 
cording to Eq.© with Ti nt 3> l/cr leads to a Kronecker 
delta 8 uu i and therefore makes any double-sum over fre- 
quencies into a single-sum. In the large- volume limit and 
for d = 1, we then obtain 



(I)(t) 



1 



2Ti n i Jq 



dio2iT (w|/ s (w)| 2 (l + coswt) + 

uj\f lo (Lo)\ 2 (l-costoT)) . (20) 



If Cd s ,Cdi 3> then, within a good numerical approx- 
imation, we can first replace the linear to dependence 
in Eq. (|20[) with Co s and Cdi in front of the correspond- 
ing angular-frequency distributions, and then extend 
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the integration to include arbitrarily negative angular- 
frequencies. One then finds that 



W(o) 



COS TU! S 



UJlo 



■ COS TUllo 



(21) 



where the interference effects are exponentially sensitive 
to the spectral width a of the single-photon angular- 
distributions similar to the spectral width dependence 
in the famous Hong-Ou-Mandel two-photon experiment 
[27| . The asymptotic value of (7}(r)/(7)(0) is given by 
(l + w /o /w s )/2. In FigHwe exhibit (7)(t)/(/)(0) accord- 
ing to Eq. (f2"0)) with the choice as in Eq. (fT5]) . It is now 
clear that with a given reference distribution \fi (u})\ 2 
of the local oscillator one can, e.g., infer the common 
spectral width a of the single-photon sources as well as 
the corresponding angular- frequency u) s . A related ex- 
perimental situation is discussed in Ref. [3l[ for a general 
one-photon state, i.e. not necessarily a pure quantum 
state. 

In the case of coherent states in the signal and local 
ports, with one-particle state angular-frequency distribu- 
tions / s (w) and fioioj) as above, Eq. (fT9|) is now modified 
according to 

E {+) {t)\i) = iJ2 \sujuo 1/2 (if s (ujy Mt \l + e WT ) 

+/ Jo (w)e < *»<*>(l - e iUT )) \f s ) c ® \f lo ) c , (22) 

i.e. Eq. ([2H|) is, for real- valued single-photon distributions 
/ s (w) and // (w), replaced by 



(I)(r) = 



2 



dujL}2n (|/ s (w)| 2 (l + cosujt) + 



|/; (o;)| 2 (l - coswr) - 2/ s (cj)// (w) sinwr) . (23) 

In Fig(3] we exhibit (I) (r) / (I) (0) with coherent states 
generated by the same choice of single-photon states as 
in Fig 121 The interference is, as in Fig 121 sensitive to the 
actual angular-frequency distributions. With a given ref- 
erence distribution |/; (a;)| 2 of the local oscillator, one 
can now infer the common spectral width a of the coher- 
ent state sources as well as the corresponding angular- 
frequency ui s . By a comparison with Figj2l we conclude 
that the MMI setup is sensitive to the actual form of the 
initial quantum states despite the fact that we are only 
considering single-photon detection processes. 



IV. INTERFERENCE OF THERMAL LIGHT IN 
THE MICHELSON-MORLEY 
INTERFEROMETER 

A. Thermal light 

For the readers convenience, let us first outline a sim- 
ple and quantum mechanical description of thermal black 




rcr 



FIG. 3: The normalized intensity (7)(r)/(/)(0) as a function 
of the time-delay r/cr for the case of coherent states as gener- 
ated by the same one-photon states as in Fig[2] Two different 
examples are plotted with 0j s = 3a: uji — 3.15a with (solid 
line) or uji = 2.85<r (dashed line). The asymptotic values of 
(-0 M / (I) (0) are as in FigH 



body radiation at an absolute temperature T. Black 
body radiation is, basically, light with random phases. 
For a single mode with angular-frequency u) and if \nj) = 
a? (u)) n " \0) / y/n u \ denotes a n w -photon state, the density 
matrix describing the thermal light is given by 



in terms of the Bose-Einstein distribution 
n(uj,T) Y 1 



Pn(w = , =rl 

V 1 + n(uj, 1 )J 



with 



n(w,T) 



+ n(u,T)J l + n(w,T) 



1 

exp(Tkj/kBT) — 1 



(24) 



(25) 



The state p{oS) corresponds to an extreme value of the 
von Neumann entropy S in units of fcgT, i.e. 

S=-Tr[pHlnpH], (26) 

subjected to the constraints (see e.g. Refs.flll [l9|) 



Tr[pH] = 1, 
Tr[a'(u>)a(u})p(u))] = n(iv,T) . 



(27) 



The random, or chaotic, nature of the quantum state 
p(u) corresponds to a phase-independent Glauber- 
Sudarshan 'P(a)-representation (see e.g. Refs.[32l-[35|) in 
terms of a coherent state, i.e. 



p(u) = / d 2 aV(a)\a)(a\ , 



(28) 
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using a single- mode coherent state \a) = exp(aa^(a;) 
a*a(u)))\0). For thermal light one finds that 



P(a) 



1 



■exp(-\a\ 2 /n(uj,T)) 



7m (cj, T) 

which obeys the normalization condition 
Tr[p(ui)] = J d 2 aV(a) = 1 , 



(29) 



(30) 



as well as 

Ti[a}(u])a(uj)p(u>)] = I d 2 aV(a)\a\ 2 = n(ui, T) . (31) 



A multi-mode system at thermal equilibrium is then 
described in terms of a tensor product p(T) = p(u>), 
where we have performed the replacement a — > a(u>) 
in Eqs.|J2SJ| and ((29}. The Glauber-Sudarshan V{a)- 
representation for the state p(T) is now, in particular, 
useful in our considerations since the response in a single- 
photon detector can be obtained immediately from the 
previous results for coherent light using an average pro- 
cedure. 



Particular values are J(l) = ir 2 /6 and J(3) = 7r 4 /15. 
In the case of d = 3 we, therefore, recover the well- 
known expression for (7)(r)/(7)(0) [19| as also discussed 
in Ref.[l_8j. In, e.g., d = 3 it is actually possible to carry 
out the relevant integral in the expression Eq.([34|) for 
(7)(r)/(7)(0) exactly with the result 



(JM 
(I)(0) 



1 + 15 



+ cosh(2a7r) 3 
sinh(a7r) 4 ( a7r ) z 



(36) 



Eg. (136)) shows that interference effects have a power-law 
sensitivity for larger a. In Fig[J]we show (7)(t)/(7)(0) 
for varying values of a in the case with d = 3 and one 
infers a characteristic coherence length l c of thermal light 
in the MMI of the form 



L~1.5 



he 



(37) 



as also discussed in Rcf. 181. 



V. THERMAL LIGHT IN THE SIGNAL AND 
LOCAL PORTS 



B. Thermal light in the signal port and vacuum in 
the local oscillator port 

We are now in the position to consider a density matrix 
describing a system where we have thermal light in the 
signal port and vacuum in the LO port, i.e. the initial 
density matrix p of the total system is given by 



p = p(T)®(\0)(0\) lo . 



(32) 

By making use of Eg. (j2"3")) . with f s {uj) — > a(ui) and 
fio(u}) — > 0, and then performing an average over a{uS) 
according to Eq. (|31[) . we immediately obtain the result 



(I)(r) 



1 



2T, 



int JQ 



duj2iruj d n{uj,T)(l + cosLUT), (33) 



in d space dimensions. In passing, we notice that 
Eg. ([33]) is actually valid for any physical quantum state 
of the form Eq. (13"2")) due to the generality of the Glauber- 
Sudarshan 7 , (a)-representation Eg. ([2"B|) . The interfer- 
ence effects as exhibited by the MMI setup therefore 
only depends on the, in general, angular-frequency de- 
pendent mean-number Tr[a^ (ui)a(uj)p\ — n(u) and not 
on other features of the actual quantum state considered. 
By a straightforward change of integration variables in 
Eg.p3p. and with a = rksT/h, we then find that 



1 



1 



J(d) 



' / x d cos(ax) 

dx \ — r\ — t 

Vexp(a;) — 1 



(34) 



where J(d) can be expressed in terms of gamma and Rie- 
mann's £ functions, i.e. 



j[d) = r(i + d)((d + 1) = r(i + d) ■ 



i 



,d+i 



(35) 



As we have seen above, with a vacuum in the LO port 
and with a signal thermal source the single-photon de- 
tection process exhibits an interference pattern. We now 
investigate what happens if we have thermal light with 
a temperature T in the LO port and thermal light with 
temperature T\ in the signal port. For this setup the 
corresponding initial density matrix becomes 



p = p{T x ) ® p(T ) . 



(38) 



By making use of this density matrix as well as the same 
methods as described in Section IIVBI by performing in- 
dependent averages over f s (tu) and /z D (w) in Eq. (|23[) ac- 
cording to Eg. pip , we immediately obtain, in the large- 
volume limit, the result 



(I)(r) 



1 



2T, 



int Jo 



dujui d 2ir (n(w,Ti)(l + coswr) + 

n(w,T )(l -coswt)) . (39) 



Here we notice that the last term in Eq. (|2"3"]) . suitably 
extended to complex- valued and // Q (w), averages to 

zero due to the chaotic nature, i.e. phase-independence, 
of thermal light according to Eqs. ([28]l and (p9|l. The 
relative intensity with respect to zero time-delay r has 
then the form 



(i)(0) 



fTps4 



15 

2^4 



dx 



exp(x) 



— (cos(aiir) - (p^) 4 cos(a a;)) 



(40) 



where a = TkgT /h and a\ — rkBTi/h. The integrals 
in Eo. (|4U| can, again, be solved analytically in a fashion 
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FIG. 4: The normalized intensity (7)(r)/(/}(0) as a function 
of a = rkBTo/h for the case of thermal light in the signal 
port and vacuum in the LO port and for d — 3. We infer a 
characteristic thermal coherence time r c in terms of a ~ 1.5 
(dashed vertical line), i.e. r c ~ l.^h/ksT. 
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FIG. 5: The same as in Fig|4]but with thermal light in both 
ports with Xi = 1.01Tb (dashed-dotted line), Ti = 0.99T 
(dashed line), and Ti = To (solid line). The vertical dashed 
line corresponds to the thermal coherence length l c = ct c with 
TckgTi/h ~ 1.5. 



similar to the integral in Eq. dMl) . i.e. 



(I)(r) 1 
(7>(0) 2 



/?0n4/ 2 + cosh(2a 7r) - 

^Ti' sinh(a n 7r) 4 • 

15 2 + cosh(2ai7r) 



2 sinh(ai7r) 



(41) 



Due to the presence of hyperbolic functions in Eq. (|41~|) . 
we observe that ( I) (r)/(I) (0) approaches its asymptotic 
value (1 + (Xb/Ti) 4 )/2 exponentially fast as a function of 
the time-delay r in contrast to the power-law dependence 
in Eq. ([36| . It is now of interest to study the behavior of 
Eq. (l4"TT) when the temperature To and T\ of the local 
oscillator and signal respectively are varied. In Figj5]we 
exhibit the interference when To is slightly smaller or 
larger than T\ as a function of the parameter rksTi/h 
for d = 3. The corresponding interference, of course, 
disappears when the two temperatures are equal. The 
sensitivity of the interference pattern with regard to the 
difference in temperatures of the source and the reference 
temperature, i.e. of the local source, constitutes the basic 
ingredient of the FIRAS setup. We also observe that the 
coherence length l c for this MMI setup is the roughly 
same as in Section HVBl 



VI. CONCLUSIONS AND FINAL REMARKS 

We have seen that the interference of thermal light in 
the Michelson-Morley interferometer can be described, 
in a straightforward manner, by making use of Glaubers 
theory of photon detection and elementary quantum 



theory of the electro-magnetic field. Furthermore, we 
have seen the emergence of a natural coherence length 
l c ~ Tic/ksT of thermal light in the MMI. It may be sur- 
prising that non-monochromatic and chaotic light, with 
random phases, exhibits interference effects but, as we 
have seen, such an interference is naturally traced out 
in terms of the normal-mode expansion of the quantized 
electro-magnetic field (see e.g. the comments in Ref.|36]). 
The result for thermal light in both the signal and the 
local oscillator ports shows that the interference pattern 
is sensitive to the difference in temperature of the two 
sources. This is the basic principle used by the FIRAS on 
board the COBE satellite in order to perform high preci- 
sion measurements of the temperature and the spectrum 
of the cosmic micro-wave background radiation. 

Since we have been considering initial quantum states 
in terms of a fixed number of photons as well as "classi- 
cal" states, corresponding to coherent states with an in- 
finite number of photons present, a quantum-mechanical 
language is mandatory. The signal and local ports of the 
MMI setup corresponds to independent input sources. 
It is, of course, a well-known experimental fact that in- 
dependent photon sources can give rise to interference 
effects (see e.g. Refs.[13-S2|)- Despite the fact that 
such interference effects are well established, the inter- 
pretation of them can, nevertheless, give rise to interest- 
ing issues regarding the very fundamental aspects of the 
quantum- mechanical world (see e.g. (43l |) when consider- 
ing, in particular, interference effects using single-photon 
sources. 

We have seen that for multi-mode systems the quan- 
tum nature of theses independent sources actually effects 
the nature of the single-photon intensity measurements. 
We have already mentioned that the angular-frequency 
distribution of a single photon can be measured using a 
similar experimental setup as the MMI considered in the 
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present paper [3l|. With a vacuum state in the local os- 
cillator port and a signal single-photon angular-frequency 
distribution /(to) of the form considered in Eq.(j9|), one 
finds, using Eq. (|2"0)) , that 

iH( i+eH " >VWo -)' (42) 

We conclude, in view of our considerations, by notic- 
ing the characteristic exponential behavior in Eq. (|4"2"|) for 
single-photon interference is not necessarily the indica- 
tion of a Lorentzian angular-frequency distribution as 
claimed in the literature when, e.g., measuring the photo- 
luminenscene signal of a single quantum dot [44J or in 
studies of other interference effects of dissimilar photon 



sources (45j | . 
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